function []=display(s)
  disp('value:');
  disp(s.v);

  disp('deriv:');
  nd = size(s.d);
  if(nd(2)==1)
    disp(reshape(s.d,[nd(1) nd(3:end)]));
  else
    disp(squeeze(s.d));
  end
